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We discuss the emission of high momentum lepton pairs in central Au+Au collisions at RHIC 
(Viivjv =200 GeV) and Pb+Pb collisions at the LHC (^/sITn =5500 GeV). Yields of dileptons 
produced through interactions of jets with thermal partons have been calculated, with next-to- 
leading order corrections through hard thermal loop (HTL) resummation. They are compared to 
thermal dilepton emission and the Drell-Yan process. A complete leading order treatment of jet 
energy loss has been included. Jet-plasma interactions are found to dominate over thermal dilepton 
\Q ' emission for all values of the invariant mass M. Drell-Yan is the dominant source of high momentum 

lepton pairs for M > 3 GeV at RHIC, after the background from heavy quark decays is subtracted. 
At LHC, the range M < 7 GeV is dominated by jet-plasma interactions. Effects from jet energy loss 
on jet-plasma interactions turn out to be weak, but non-negligible, reducing the yield of low-mass 
» I ' dileptons by about 30%. 



43 



43 



X 

5-H 



I. INTRODUCTION 



< 
(N . 

Finding experimental evidence for the existence of a Quark-Gluon Plasma (QGP) is one of the main reasons for 
£N) ' conducting relativistic heavy ion collisions. In this context, real and virtual photons can be used to probe the high 
temperature and density phase in these collisions: due to their large mean free path, electromagnetic radiation suffers 
only little final state interactions Q, 0. More recently, photons emitted from jets interacting with the medium have 
been suggested as a further signature^ . A quantitative characterization of the different sources contributing to the 
photon yield in high energy heavy ion collisions can be found in [4j . In this work, we concentrate on virtual photons 
decaying into lepton pairs. 

Several sources for lepton pairs compete and have to be considered: dileptons from the Drell-Yan process thermal 
dileptons [|| both from the QGP 0, M, and from the subsequent hadronic phase 0, 0] , dileptons from the absorption 
of jets by the plasma and dileptons from bremsstrahlung of jets. Besides the thermal emission, dileptons from 
the latter two sources carry information about the medium. Comparing with the baseline p + p process, we note that 
absorption of jets exclusively happens in nucleus-nucleus collisions in the presence of a medium. There is a also a 
significant change in bremsstrahlung emission expected in a medium (in A+A) compared to the case of vacuum (in 

p + p ">- — . 

The main background for dileptons are correlated charm decays (DD — » e + e~ X) at intermediate mass 12], while 

one has to cope with Dalitz decays of light mesons (ir a — > je + e~ , u> — > 7r°e + e~) at low mass M < 1 GeV. In principle, 

there is also another source of dileptons corresponding to pre-equilibrium emission. It is difficult to assess both 

theoretically and experimentally, nevertheless, such calculations have been attempted ^| . 

In Ref. the dilepton yield from the passage of jets through a plasma has been evaluated as a function of invariant 

mass at leading order. The results show that this process dominates over thermally induced reactions at high invariant 

mass M . However, energy loss of jets in the medium has not been included in this calculation. Large energy loss 

of jets was discovered as one of the most exciting results from the Relativistic Heavy Ion Collider (RHIC) and has 

been observed in the suppression of high px hadron spectra |l4l Il5| and through the disappearance of back-to-back 

correlations of high px hadrons [0| . 

In this paper we re-examine the leading order dilepton production from jets an d we explore the effect of energy 

loss. We assume here that jets lose their energy by induced gluon bremsstrahlung only 0,0. We use the approach 

developed by Arnold, Moore, and Yaffe (AMY) |20|. which correctly treats the Landau-Pomeranchuk-Migdal (LPM) 

effect (up to 0(g s ) corrections) j^]. This model has been used successfully to reproduce the measured nuclear 

modification factor Raa of neutral pions, the ratio of all photons over background photons and the direct photon 

yield at RHIC Jets will be defined by all partons produced initially with transverse momentum ^> 1 GeV. 

The total dilepton production could be influenced by the choice of the cutoff p 3 ^ ■ As we will discuss below, in order 

to avoid such sensitivity, we limit our study to high momentum dileptons. 

In perturbation theory at high temperature, it is important to distinguish between hard momenta, on the order 

of the temperature T, and soft momenta, on the order of g s T, where the QCD coupling constant is assumed to be 

small, 3 S < 1. Only at the end of the calculations, will the results be extrapolated to more realistic value of the 

coupling constant. When a line entering a vertex is soft, there are an infinite number of diagrams with loop corrections 



contributing to the same order in the coupling constant as the tree amplitude. These corrections can be treated with 
the hard thermal loop (HTL) resummation technique developed in Ref. [22 . This technique has been used thereafter 
in Refs. jH, 0> |2]| to show that the production rate of low mass dileptons (M < T) could differ from the Born-term 
quark-antiquark annihilation by orders of magnitude. In this work, we apply this technique to go beyond the leading 
order jet-medium interaction. 

The dilepton production rate in finite-temperature field theory is calculated in Sec. |H] and the physical processes 
underlying each contribution are discussed in Sec. IIIII in the framework of relativistic kinetic theory. In Sec. IIVI we 
include these dilepton production rates in a model for a longitudinally expanding QGP fireball. In Sec. we calculate 
the yield from the Drell-Yan process and correlated leptons from heavy quark decays. The numerical results are 
presented in Sec. I VII Finally, Sec. IVIll contains our summary and conclusions. 



II. DILEPTON PRODUCTION RATE FROM FINITE-TEMPERATURE FIELD THEORY 

Within the thermal field theory framework, the production rate of a lepton pair with momenta p + and p- is given 
by HJ 

where is the finite-temperature retarded photon self-energy. After integration over the angular distribution of 
the lepton pair, we get 

dR e+e ~ _ 2F dR e+e ' 2a dR^' a Imnff" 

d 4 p dM 2 d 3 p 3irM 2 d 3 p 12tt 4 M 2 1 - e E / T ( ' 

where p = p" + +p~L,E and M are the momentum, energy and invariant mass of the virtual photon. 

In the HTL formalism, the leading order resummed photon self-energy diagram is shown in the left hand side of 
Fig. ^ The heavy dots indicate a resummed propagator or vertex. The second diagram in Fig. ^ coming from the 
effective two-photon-two quark vertex, is needed in order to fulfill the Ward identity 

(p) = o. (3) 

Using power counting [22L |27| . the HTL resummation gives corrections of order <? 2 T 2 /|p| 2 to each bare vertex of the 
first diagram of Fig. ^ while the bare propagators receive g 2 T 2 /\q — p\ 2 and g 2 T 2 /\q] 2 corrections. A resummed 
propagator consists of an infinite number of gluon correction, as shown in Fig. [21 

In this work, we are interested in the high momentum dilepton limit 3> T), so that vertex corrections can be 
neglected, and at least one propagator is guaranteed to be hard. We assume that this is the propagator associated with 
momentum q — p. Then the propagator associated with q can be soft or hard, implying that is has to be resummed. 
Also, for the high |p| limit, the second diagram of Fig. Ogives only <? 2 T 2 /|p| 2 corrections to the first diagram. The 
resulting diagram that has to be evaluated is shown in Fig. [3] and its contribution to the self-energy is 

n»=3e 2 ]T(^) 2 T£ J ^TT[rS D (qh,S(q~p)} (4) 

with the sum running over the Matsubara frequencies. The diagram in Fig. [3] includes the leading order effect, and 
some next-to leading order corrections in g s . In order to have a complete next-to leading order production rate in the 
region \p \ 3> T, contributions like bremsstrahlung need to be included. This will be discussed in the last section. 
The dressed fermion propagator, in Euclidean space (gf = —iqo"f° + q • 7), is given by [2^| 

Q f s 1 7°-7-q , 7° + 7-q fK ^ 

SDiq} = 7T^ = ^qy + ^irw (5) 

where 

D±(q) = -iq ±\t\+A±B (6) 
The terms A and B describe the quark self-energy 

S = V + Bj ■ q. (7) 
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FIG. 1: Resummed photon self-energy diagrams. 
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FIG. 2: Resummed quark propagator. 



In the HTL approximation, one obtains [27 



(8) 



where mp — gsT/V6 is the effective quark mass induced by the thermal medium and the Q n are Legendre functions 
of the second kind. The bare propagator 5(g) is the same as 5u(g), with D± replaced by 



d ±(q) = -i<io ± |g1 

Using those expressions for the quark propagators and carrying out the trace, Eq. pleads to 

d 3 q 

1 /l-q-k l + q-k\ 1 / 1 + q ■ k l-q-k" 



(9) 



D+{q) \ d+(k) d-(k) / D-(q) \ d+(k) d-(k) 



(10) 



where we have defined k = q — p. 

At this point it is convenient to introduce the spectral representation of the effective quark propagator defined by 



P±(u, \q\) = -21m 
1 



1 



D±(iq , |g|)' 
dw p±(w,\q\) 
D±{iqo,\q\) J 2tt iq - u + it ' 



This implies 



P±(w, |gl) 



2tt 



2m% 



[S (w - w±(|$|)) + <5 (w + c T (|gl))] - 2tt/3 ± (u, \q\) 6 ( |gf - 



(11) 
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FIG. 3: Effective photon self-energy diagram for hard external momentum p. 
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with 



/3±(^|gl) = 



-u±\q\ 



\<T\ 



(±1 + | In 



+71 



-I '71 



(13) 



Here w± = w±(|<z|) correspond to the poles of the effective quark propagator. They are solution of D±(u>, \q\) = 0. 
The solution to = oj + represents an ordinary quark with an effective thermal mass y/2m p and a positive helicity over 
chirality ratio, x = 1 |23(. The solution ui — lu- represents a particle having negative helicity over chirality ratio, 
X = — 1- This collective mode, called plasmino, has no analog at zero temperature. Following the notation of Ref. [23| . 
we denote the ordinary modes by q+ and the plasminos by q_ . The spectral density of the bare propagator is simply 



r±(uj',\k\) = -21m 



1 



d±(ik , \k\) 



-2ttS(lo't \k\) 



(14) 



Using the spectral density representation has the advantage that one can carry out the sum over Matsubara frequencies 
with help of the elegant identity |23( , 



ImT 22 Fx (iq )F 2 (iq - ip ) 



IT 1 



e E/T) I 

, 27T 



dtu' 

~2r7 



xpi(aj)p 2 (-uj')S(E -u>- cj')/fd(w)/fd(w') 



(15) 



where /fd is the Fermi-Dirac phase space distribution function and pi is the spectral density associated with Fi . We 
use the analytical continuation ipo — > E + it with the dilepton energy E. 
Putting all information together, we obtain 

tan »_^j:(?)'(i-^)/^£|A»M 

.(w, \P\){S{E - u + £ x )(l - q- k)f FB (-Ei) +S(E-cu-E 1 ){l + q- k)/ FD (ia)} 
.(w, \P\){S(E -U3 + S x )(l + q • k)f FD (-Ei) +S{E-u-E 1 )(l-q- k)/ FD (£i)} 



-P- 



(16) 



with Ei = \q — p\. Since both E ^> T and E\ 3> T, the terms proportional to <5(i? — u> + E\) are exponentially 
suppressed by the factor / F d( w ) for u> 3> T. Also, because £i corresponds to a massless parton, energy and momentum 
conservation does not permit the terms containing 6(E + \ui\ — E\) for lo 2 > \q\ , i.e. there is no phase space available 
for the process 1 — > 2 + 3 if parton 1 is on-shell. These arguments lead to further simplifications and we have 



imir?" 



d 3 q 
(27T) 5 



/fd(w+)/fd(£'i)- 



2m| 



■ ( y(S-w+-£; 1 )(l + q-k) 



+/fd(w-)/fd(-^i)- 



2mf, 



<5(£-ct>_ -£i)(l -q-k) 



■ / da;/ FD (a;)/ FD (S 1 ){/i + (a;J^)(l+q-k)+^_(a;,|^)(l-q-k)}5(£-a;- J B 1 



(17) 



Now we analyze the integral over d 3 q = \q\ 2 d\q\ dQ,. When \q\ becomes as large as E\, the propagator associated 
with q does not have to be resummed. If \q\ gets much bigger, implying E\ ~ g s T, then this is the propagator 
associated with E\ which should be resummed. We have a symmetry relatively to the point \q\ = E\, such that 
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FIG. 4: Processes extracted from photon self-energy in Fig. [3] In the picture, the particle with energy E\ is an anti-quark, 
but the quark contribution is also included in the calculation. 



Jo°° ^ 1^1 ~ * \T ^ I Si 2©(-Ei — |<z|)- According to Eq. @ the dilepton production rate is thus given by 



dR e+e 
d 4 P 



(27T) 



/fd(w+)/fd(£'i) 



<4 - \Q\ 2 



2m| 



<5(£ - w+ - £?i)(l - q« pi)0(£?i - kl) 



+/pd(w_)/fd(-Ei)- 



2m| 



<5(£ - u_ - £7i)(l + q - pi)e(£7i - 



+ / dufp D (w)f PD {E 1 )5(E -u- Ei)Q(E! - 

— OO 

: {jS+Cw, |$|)(1 - q ■ Pi) + /9_(w, + q • Pi)} 



(18) 



where pi = pi/|pi| = — k. 

The dilepton production rate from finite-temperature field theory calculated in this section follows the path pre- 
sented in Refs. jl^,|2^|. New to our approach is that we keep the general form of the function /fd(Ei) until the end of 
the calculation, while this function is integrated out and information about it is finally lost, in the earlier work. In the 
next section we will analyze, for the first time, the physical processes behind the production rate in relativistic kinetic 
theory. This is done so that /fd(-E'i) can be interpreted as the phase-space distribution function of an incoming parton 
of energy E\ . Once this is established one can apply the usual technique and obtain the production rate involving 
jets by substituting the thermal distribution /fd(-E'i) by the time-dependent jet distribution /j et (-Ea, i). 



III. DILEPTON PRODUCTION RATE FROM RELATIVISTIC KINETIC THEORY 

We will show in this section that the production rate in Eq. (|18l) can also be interpreted in relativistic kinetic theory. 
Indeed, after cutting the diagram in Fig. [31 we get the Feynman amplitudes shown in Fig.^J The process in Fig.QJa) 
corresponds to the annihilation of the hard antiquark of energy E\ with a soft quasiparticle q+ or q_. This is the 
pole-pole contribution, as it involves the propagation of the poles of the propagator in the diagram shown in Fig. [3] 
Fig. 0{b) corresponds to Compton scattering of an antiquark with energy E\ with a hard gluon from the medium, 
with the exchange of a soft quasiparticle and Fig. 01 c) represents the annihilation of an antiquark E\ with a hard 
quark from the medium with the exchange of a soft quasiparticle. They correspond to a cut through the self-energy of 
the dressed propagator in Fig. [31 this is thus called the cut-pole contribution. There is no s-channel process because 
hard thermal loops have imaginary parts for q 2 < only .27]. The dots mean that the quasiparticle propagator is 
resummed. That modification of the propagator in the space-like region at non-zero temperature is known as Landau 
damping. 

The production rate, as computed from relativistic kinetic theory (kt), for reaction 1 + 2— > 7* + 3 + ... is 

4=/(n^fe) <* w» — -•) ^ ' ( *»" a.) 

where M. is the invariant scattering matrix element. 
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A. Pole-Pole Contributions (1 + q± — ► 7*) 

The production rate for the process of Fig. Ufa) is 



E 



kta 

d 3 p 



d 3 Pi d 3 q 4 4 . .2 /*+«~(£i)/fd(<7o) 

(2tt) (5 (pi + g-p) |M kta | - 



2(2tt) 3 £i 2(27r) 3 g 



2(2tt) s 



8(27r)5g £] 



5(E 1 +qO-E) \M kt J / 9+9 (£i)/fd(<7o 



(20) 



where f q+q (E\ ) is the phase-space distribution of quarks and antiquarks with energy E\ . The square of the matrix 
element, after summation over spin and color is 



M 



kta I 



\Mi 



6e 2 



Tr 



Mi +q _-, r * I 

#1 (£ < + (?)< + (?) + E <~ (9) fi 2 _ (?) 



(21) 



We need to find a completeness relation for J2 s (g , )ul ± (?)• The dressed propagator in Minkowski space is 

1 



and by definition, it is |3flj 



iS D (q) 



iS D (q) = 



7° - 7 • q 7 + 7 • q 



rf-E 2#+(<z) 2Z?_(<z) 

E s ^ + (<z)^ + (<z) , E.«2 _ (9)fi2-(9) 



<? 2 - 



<? 2 — m 2 _ 



(22) 



(23) 



where m+ and m_ are the physical mass of the quasiparticles. We expand D± around the pole located at qo = u>±: 



D±(q) ~ {q - uj± 



dD± 



= («b-w±) ( -l + ^-(A±B) 
oqa 



m 2 F {q 2 - m|) 



From Eqs. (|2*2*j) . l|2*3)l and we finally obtain the relation 

2^ M * (?M (?) = 



2m 2 , 



(7 Tr q) 



Substituting this into Eq. (|21|l gives 



\M 



12e 2 2 (f£) 2 4- [ w+ (^ 2 - | 9 f)(l - q • Pl ) + cj_(oj 2 _ - \q\ 2 )(l + q • fc)] 



kta I 



(24) 



(25) 



(26) 



and the virtual photon production rate becomes 



di?7, 12 
1 & 



d 3 p 



111 



d-- c 



8(2tt) 5 

{fM"+)f 9+ HEl)(ul - \qf )(1 - q • Pl )5(i?l + c+ - £)6(^ - 

2x 



+/ FD (u;_)f' + «(£ 1 )(u; 2 _ - |gf )(1 + q ■ p 1 )5{E 1 + cu- - E)0(E 1 - \g\)} (27) 
We have introduced here the same cut Q(Ei — \q\) as in the previous section. Finally the dilepton production rate 



G 



from the process in Fig. Ufa) is 



dRf + , e 2a dRf, 

ktn _ AUL tp Kta 

d 4 p 3nM 2 d 3 p 

2 2 r°° i^i 2 



tt^^UJ y o (2tt)3 y 

x {/ FD (^ + )/^(^ 1 ) ( ^+ TO j g1 } (1 - q • p 1 )S(E 1 +uj + - E)0(E 1 - \q\) 

+fMu-)f q+ HE 1 ) {UJ - 2 Jf \ l + q • p^E, + uj_ - E)Q(E 1 - \q\)} . (28) 



B. Cut-Pole Contributions (1 + 2^3 + 7*) 

The expressions for the annihilation and Compton scattering processes from Figs. 0Jd and 0J; are 

J R kUc f d* Pl d 3 p 2 d* P3 | 2 /^(£ 1 )/be(£2)(1-/fd(£ 3 )) 
E ^p~ J 2(2^2(2^ 2(2nfE 3 {2n) 6 ^+P*-P-P*)\ M kt b \ ^3 

d 3 Pl d 3 p 3 d 3 p 2 4 4 , |2/^(S 1 )/ FD (^ 3 )(1+/be( j B2)) 

-(2tt) d (pi +p 3 -p -p 2 ) |A1 ktc ' 



2(2^)3^ 2(27r)3£;3 2(2tt)3£; 2 v ; > r* r r*J | — Ktcl 2 ( 2?r )3 

(29) 

where /be is the Bose-Einstein distribution function. The squared matrix elements of the diagrams in Fig. and 
Fig. are given by 

i^kt 6 r=i^kt c r= 3 ^ 2e2 E(-) 2Tr [ 4 ^^^-^)^^b-pi)] , (30) 

where Cf = 4/3 is the quark Casimir. It is important to point out that only the t-channel exchange is considered 
in the Compton scattering process. For the annihilation process shown in Fig. the particle with energy E\ is an 
antiquark and the one with energy E 3 is a quark. There is also a contribution with E\ and E 3 being associated with 
a quark and an antiquark respectively. Those two contributions are added incoherently since coherence effect are 
suppressed by higher powers in g s . After inserting 1 = J_ du> S(ui — E + E\) and using p\ = p — q we obtain 

E dR ^ r du: f ^ ^ 3(x E + F ) 

S(uj -E 2 + E 3 ) \M ktb \ 2 f q+ ^E{)f B ^E 2 )[l - /fd(£ 3 )) 
+5(oj -E 3 + E 2 ) \M Vic \ 2 P+HE^hniE^l + f BE (E 2 ))} 

To compare this result with the one obtained in the last section we take advantage of the quark self-energy E, 
shown in Fig. |3J The expression for the discontinuity of E in the space-like region is | .3l| 

Disc E(w, \g\) = -mC F g 2 s J 7^4^ - E 3 + £ 2 )/fd(£ 3 )(1 + /be(£ 2 )) 

+6(lj + E 3 - E 2 )f BE (E 2 )(l - /fd(£ 3 ))}/fd 

(ui) (31) 

Eqs. (EOJ, (EU and ^ lead to 

x Tr [4 faS* D {q) Disc E(w, \q\) S D (q)] (32) 
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FIG. 5: Quark self-energy with gluon propagator in the loop. 

We can use the relation 

^,(g)DiscE(w,|5l)5 D (g) = Disc(-<5 D ( g )) (33) 

which holds given that D±(q) = D±(q*). This is indeed the case as can be inferred from the definition of D±, Eqs. © 
and ijSJ. For uj 2 — \q\ 2 < 0, we use Eqs. (fTTJl . (|T2*|l and j^H to express the right hand side as 

Disc (-<&,(<,)) = - ^-J-^ Disc^- ^+J-^ Disc 1 



2 ^ £>+(g) 2 ~ £>_(<?) 

-1(7° — 7 ■ q)Im — — i(jy + 7 • q)Ir 



D+(q) w ^ £>_(g) 

= -z^(7 - 7 .q) / 3 + (a;,|gl)- i 7r(7 +7.q) / 3_( W ,|< 7 l). (34) 

Using the latter result in Eq. I|32l) and carrying out the trace, we find that 

X [P + (tu,\q^)(l - q - p\) + p^cu,\q\)(l + q- p 1 )]0(E 1 ~ \q\) (35) 

As before, we have introduced the term Q(E\ — \q\) as we consider only the region where HTL may be important. 
The dilepton pair production rate for the process shown in Figs. 0J) and c is 



d 4 p 3ttM 2 d 3 p 



2a2 ^(^) 2 f du, P^Jd^-E+E^mfM-) 



n 2 M 2 



x [/3+(w, \q\)(l - q • pa) + p-(w, \q\)(l + q • p\)] e(E 1 - \q\) (36) 

Upon adding Eqs. (|2*%|) and we reproduce the result from Eq. (fH^I . when the particle associated to E\ is thermal, 
i.e f q+ ^(Ei) — ► 2/fd(^i)- This proves that both methods, finite-tempcraturc field theory and the relativistic kinetic 
formalism, lead to the same result. 

We now briefly compare our approach with the method used by Thoma and Traxler in Ref. |25|. They have 
calculated the photon self-energy shown in Fig. [21 with an imposed cutoff k s <C T on the momentum \q\ in the loop- 
integral, such that < \q\ < k s . They then added the Compton scattering and annihilation processes coming from 
cutting the two-loop photon self-energy without HTL propagators or HTL vertices. Those two latter process have an 
infrared divergence, which is regulated by imposing a low value cutoff k s for the exchange momentum. When adding 
all those processes, the final production rate is infrared safe and independent of k s . They have also calculated the 
a 2 a s contribution coming from the pole of the effective quark propagator in Fig. [31 However in their approach, the 
information about the parton phase space distribution is lost, i.e. it is not possible at the end to make the substitution 
fq+q _> Here, we only consider the one-loop diagram from Fig. [21 but we use the dressed propagator So(q) up 

to the scale \q\ = k c , where k c corresponds to E\ due to the 9(Ei — \q\) function. With this method we do not have 
to specify the shape of / 9+<7 until the end of the calculation. We have verified that our numerical result depends only 
weakly on the scale k c . For example, taking k c = 0.6 x E\ reduces the production rate by ~ 20 %. 
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FIG. 6: (Color online) Production rate of dileptons with momentum p — 4 GeV, from thermally induced reactions, at a 
temperature T = 300 MeV and for a s = 0.3. Dotted line: Born term; short-dashed line: pole-pole contribution with particles 
having negative identical helicity over chirality ratio \\ dot-dashed line: pole-pole contribution with particles having identical 
helicity over chirality ratio Xi long-dashed line: cut-pole contribution; solid line: sum of all processes; and double dot-dashed 
line: Born term plus a 2 a s contributions from Ref. |2fi|. 



Fig El shows, for f q+ ^{Ei) — ► 2/ FD (£? 1 ), the different sources of dileptons at a temperature T = 300 MeV. In all 
cases, the particle with energy E\ corresponds to a pole with positive %■ The pole-pole contributions are shown by the 
dot-dashed and the short-dashed lines. They correspond to the diagram in Fig.^Ji. The annihilation of two partons 
with positive helicity over chirality ratio, \ = 1, (dotted-dashed line) dominates at high invariant mass. For M > 1 
GeV it converges toward the Born term (dotted line) obtained from a one-loop photon self-energy calculation with 
two bare propagators. The cut-pole contribution (long-dashed line) is dominant for M < 1 GeV. The corresponding 
physical processes are the annihilation and Compton processes in Fig.0Jj and c. The sum of our contribution is shown 
as the solid line. It agrees very well (within 30 %) with the sum of the Born term and the a 2 a s result from Ref. 25], 
given by the double dot-dashed line. 



IV. DILEPTON YIELD IN ULTRA- RELATIVISTIC HEAVY-ION COLLISIONS 



A. Thermal Dileptons 

High-p-r real and virtual photons are preferentially emitted early during the QGP phase, when the temperature is 
largest. Explicit hydrodynamic calculations show that the space-time geometry of the fireball smoothly evolves from 
a 1-D to a 3-D expansion [32| . By the time the system reaches the temperature corresponding to the mixed phase 
in a first-ord er p hase transition, the system is still dominated by a 1-D expansion [3^]. For such a geometry, specific 
calculations [33 suggest that the flow effect on photons and dileptons from the QGP is not large at RHIC and LHC 
for particles with transverse momentum pt > 2 GeV. 

Assuming a 1-D expansion |34| that is cut off at a maximal space-time rapidity ?7 m ax, the yield as a function of 
invariant mass M and dilepton rapidity yd is given by the rate R e e — R e e (r, 77, rj_) as 

^ = (rfirtr^VgL. (37) 



dM 2 dy d J " 'dM 2 dy d ' 

Here r, 77 and are the proper time, space-time rapidity and the transverse coordinate respectively. In the remainder 
of this work the dependence of R on r and rj_ will be omitted in the notation for brevity. As long as yd < < ?/max we 
can invoke boost invariance for dR e e /dy d to argue that 



Vn 



dR e+e ( V ) f ^ dR e+e ~(y d ) dR e ^ e ~ (y d 



drj- 



dM 2 dy 



j ''''' dM 2 dy 



v=v - dM 2 ■ (38) 
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Hence the rate of lepton pairs with a fixed rapidity yd integrated over the entire longitudinal extent of the fireball is 
the same as the rate of all lepton pairs integrated over rapidity from the slice of the fireball at rj = yd- In the following 
we study only dileptons produced at mid-rapidity yd — 0. Their yield can now be written as 



dN e+e 



dR ee On = 0) f , l2 f, f , 2tt 1 dR e e (n = 0) 
= / rdrd 2 r± >- = / TdTd 2 r ± / dp T p T / dp z — - >- . (39) 



dM 2 dy d y d =Q J dM 2 J J 1 ^ J 1 E Q 2 d 4 p 

Here Eq is the energy of the lepton pair in a local frame (where the temperature is defined) moving with rapidity 
y = |ln ^°j^° relatively to the fireball. In this local frame, the transverse and longitudinal momentum of the dilepton 

are pt and p z respectively. Then we have po — \Jp\ + p 2 and Eq = \J M 2 + p\. The dilepton energy as seen from 
the fireball is E = ^/M 2 + p 2 T . 

At this point, we add two constraints in order to facilitate comparison with experimental data. First, we introduce 
a lower cutoff PTmt for the transverse momentum of the lepton pair; second, a cut on the individual lepton rapidities 
which reflects the finite geometrical acceptance of any detector, is included, such that \y e ± \ < y cu t- We introduce a 
multiplicative factor P cut = P(|t/ e ±| < y cu t,PT) to include the latter condition. In the center of mass frame of the 
lepton pair the distribution of positive leptons, normalized to unity, is given by 



E+ dP , + - S(E+ - M) 

d 3 ptm ^Etm 



(40) 



Then the probability for a virtual photon with momentum px at midrapidity to emit two leptons with rapidities 
\lle ± I < 2/cut can be obtained by a boost back to the lab frame as 



E+ m d 3 p+ 



P(N/e±| < VcuUPT) = / "T" C? "„ 7 r " Q(l^ ± l < 2/cut) 



Note that in this frame 



7T 


\ 


E± + 


Pi 


E±~ 


P i 




M 2 



cm 

2 



d cos 0d<j) ( , 

JJ-\ 0(}y.±l<y„ t ). (41) 



M 2 

E+ = , (43) 

2( v /Af 2 +p^-pT cos 9) 

pf = ±E+sm9 sine/) (44) 



and E~ = E- E+ . 

The dilepton yield in a longitudinally expanding QGP is finally given by 



dN e e f f n± f f°° f°° 2ir 1 dR 

= / drr / drr / cbj> / dp T pr / ^ittt- P{\Ve±\ < VcuuPt) (45) 



dM 2 dy d J Jo Jo J PTmt J -oo E o 2 d 4 p 

The thermal-thermal yield at midrapidity is thus obtained by combining Eq. (|18|) and Eq. 



B. Jet-Thermal Dileptons 

In this subsection, we calculate the emission rate of dileptons from interaction of jets with the medium. The initial 
phase space distribution function for partons from jets produced in heavy ion collisions is jH l35j| 



4r(l - p -, t0) = &W£i)^ {( ,_ 9) 

9 q rpT dyd 2 p T 



(46) 



in a boost invariant Bjorken scenario. 77 is the space-time rapidity, to the formation time of the jets and r 2 =t\ — z 2 . 
dNlq I \dyd 2 pT) is the spectrum of jets and g q is the spin-color degeneracy. V{t±) represents the probability to create 
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a jet at position rj_ in the transverse plane. It is given by the product of the thickness functions of the colliding nuclei 
at rj_ normalized to 1. For central collisions and assuming hard sphere nuclei, we have 



^) = ^l-^j*(*x-r,), (47) 

where Rj_ = A x ^1.2 fm is the radius of the nucleus in the transverse plane. Since energy loss by bremsstrahlung 
involves only small 0(g s T/p) changes to the directions of particles, we suppose that jets keep propagating in straight 
lines after they are created. Then at any later time t 

41* (x, p, t) = f q f (s - , a i ) (48) 

where t = t — to is the propagation time of the jet. 

Energy loss can be described as a dependence of the parton spectrum dN 3 q ^ /dE on time. We will only be interested 
in the region around midrapidity later, so we can write 



dyd 2 p T 



= 0-1 iJL 49 

y =o EdE v ; 



with a phase space factor Q. The evolution of the energy spectrum dN%£ (t) / dE is determined by the AMY evolution 
equations. Note that we do not need to specify the constant Q to obtain the evolution of the spectrum. Here we quote 
only the basic formulas of the AMY formalism, which assumes that jets lose energy by inelastic processes only. More 
details can be found in [2(1 13(| . This formalism has been applied in Ref. 3 to successfully reproduce the ir° spectra 
at RHIC. The parton spectrum dN^/dE depends on the energy lost, but not really on how this energy has been 
lost, such that using an another model for jet-quenching, like for example a model including elastic energy loss |37| . 
should not affect our following results for the dilepton production from jet-plasma interactions. 
The evolution equations for the spectra are 



d_ f dN q f( P y 

dt 1 dE 
d fdN&ipY 



dt 



dE 



dk 



= / dk 



dMf(p+k) dr« (p+k, k) dN Q f(p) ar* (p, k) t jm et {p+k) ar?jp+k, k) 



dE 



dkdt 



dE 



dkdt 



+ 2- 



dE 



dkdt 



dN q f(p+k)dT q %(p+k,p) i dNf(p+k)dT q %{p+k,k) 



dE 



dkdt 



dNf\p) fdT^(p,k) , dT/ g (p,k) 



dE 



dkdt 



dkdt 



dE 



e(2k-p) 



dkdt 



(50) 



with the transition rates 



dT( P ,k ) _ C s gl 
dkdt 



Wirp 7 1 ± e- fe /T i ± e -(p-k)/T 



l+(l-a) 2 



(l-x) 2 
l+^ 4 + (l-x) 4 
x 3 (1— x) 3 




2h-ReF(h,p,/c). (51) 



The integration over the range k < represents absorption of thermal gluons from the QGP; the range with k > p 
represents annihilation with a parton from the QGP of energy k—p, while < k < p is the range of bremsstrahlung. 
In writing (|50|l . we used dT^ g {p,k) — dT^ g (p,p—k) and similarly for g — > qq. The 0-function in the loss term for 
g — * gg prevents the double counting of final states. In Eq. I|51|l C s is the quadratic Casimir relevant for each process 
considered, and x = k/p is the momentum fraction of the gluon (or the quark, for the case g — » qq). The factors 
1/(1 ± e~ k / T ) are Bose enhancement or Pauli blocking factors for the final states. The vector h = p x k determines 
how non-collinear the final state is, and F(h,p, k) is the solution of an integral equation describing how \p — k; k) (p\ 
evolves with time poj . 

The dileptons produced from the passage of jets through the QGP are finally obtained from Eq. (|45|l and from 
Eqs. J2gJ and (gSJ with the substitution f q+q (Ei) -> f f(Ei). Note that 77 = together with the boost invariance of 
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the jet distribution imply that the longitudinal momentum of the jet parton vanishes. After some algebra, we get 



jet— th 



dM 2 dy d 



7rM 2 g q 



/pRj_ POO pOO nOO ^ 

dt / rdr£ / dp T P{\y e ±\ < y cntl p T ) / dq T / — 
Jo J»t... Jo J-oo 



'Pt, 



2m\ 



11 W 



v/1 - cos 2 ^ ^' d yi d 2 PTl 



yi=o 



+ / «L;/fd(w)&(w, 151) 1 + i-^" - i 



|9l 



q- p 



#i M -Bi / VI - cos 2 6> d Vl d 2 PTl 



yi=o 



(52) 



with |g| = \/qx + Qz an d 



f 



4 

TTi? 2 ^ 



1 - 



r 2 + t 2 

r 2 +t 2 
Rl 



(53) 



8tr 

— - sin u , 



ttR 



for the cases that r +t — 2tr > R L , r + 1 + 2tr < R L and all other cases, respectively. Here we have defined 



wo = arccos ■ 



r 2 + t 2 - R 2 



2tr 



(54) 



The other quantities to be specified are 



q p = pt Qt cos 0± + q z , 
_i f P T + q T - (E Q - uj±) 2 



E 



2p T q T 



(55) 



E - LU± 



We obtain the jet-thermal dilepton yield without HTL effects, i.e the Born term, by keeping only the pole-pole (q + ) 
contribution in Eq. (52) • ln this case we have to substitute: 



/3±->0 : 



(u* kf ) 



2m 2 , 



o. 



4 

2m 2 , 



1 



M 2 



This leads to the final expression for the jet-thermal dilepton yield without HTL effects: 



(56) 

(57) 
(58) 



" iv jet-th 

dM 2 dy d 



^9q 



dt 



dr£ dp T j dq T / dq 



E E qy /l -cos 2 6»+ 



/ FD (£7q) j P(|7/ e ± | < y cut , pt) 



dN™(t) 



dyid 2 p Tl 



3/1=0 



(59) 



For a purely longitudinal expansion of the fireball at each point the temperature is evolving from some initial 
time Tj as 



T(r ± ,r)=T(rj. J r i )(^ 



1/3 



We assign the initial temperatures in the transverse direction according to the local density so that 0, 0, ^| 



T(r±,n) = Ti 



2 1 



R 2 



1/4 



(60) 



(61) 
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We assume a first-order phase transition and use 



/QGPM^(^-l) (62) 

as the fraction of the QGP present during the mixed phase [34[. Here r d = gc^/gn is the ratio of the degrees of 
freedom in the two phases: gq = 42.25 for a QGP with three flavors of quarks and gn = 3 for a simple gas of pions. 
Tf is the time when the temperature reaches the critical temperature of 160 MeV (see Eas. 1601 and Ifiljl . marking the 
beginning of the mixed phase, and th = r dJj is the time it ends, determined by the condition /qgp = 0. However, 
since signals associated with jets are sensitive to early times, the order of the phase transition is not crucial. The 
r-integration in is carried out from Ti to th and in addition it is scaled between r/ and th to account for the 
fact that only a fraction of the system is still a QGP, such that |j] 

T f TH 



dT= dT+ dr/ QG p(r). (63) 



V. DRELL-YAN AND HEAVY FLAVOUR DECAY 



We calculate the Drell-Yan process to order 0(a s ) in the strong coupling which is the leading order result with non- 
vanishing transverse momentum px of the lepton pair. We also have to take into account virtual photon bremsstrahlung 
from jets. The total Drell-Yan yield is the sum of the direct and the Bremsstrahlung contributions, <tdy = fdirect + 
afragllEi. 

The direct contribution in collisions of two nuclei A and B is given by 



^direct Q 2 a s \ ^ f dx a j. f _ ^ £ t _ ^ |M 0+ &_> c+7 . 



,2 



—fa/A(x a ,Q)fb/B(xb,Q) ! = r-}-= = . (64) 



dM 2 dy d dp 2 T 3M 2 s NN J x a x b a/ ' s N NX a ~ ^/s^Ny/^T p T t 

The squared scattering amplitudes |M a +b^ c + 7 * | for the Compton and annihilation processes of two incoming partons 
can be found in [io|. When px and M are both large and of the same order, the direct contribution is the dominant 
mechanism. However, when Aqcd "C M <C pt, logarithmic corrections with powers of \n(p T / M 2 ) are large. They can 
be effectively resummed into a virtual-photon fragmentation function Z? 7 » / c (z, Qf)i giving rise to the fragmentation 
contribution 



dcr frag a f dz da A+B ^ c+d 



dM 2 dy d dp\ 3irM 2 J z 2 dy d dp 2 cT 



PcT—Pt/^ 



■xD T/c (z : Q F ) . (65) 

The cross sections for the production of a massless partons c in A + B collisions, da A+B ^ c+d / dy d dp 2 cTl can be found 
in 0] . The factorization scale Q and the fragmentation scale Q p are both set to the order of the energy exchanged 
in the reaction \J M 2 + p\ 0|. The effect of varying the scale Q = k\Jp\ + A/ 2 , from fc = l/2tofc = 2, introduces 
a variation of the dilepton yield by ~ ± 35 % relatively to k = 1 . 

The typical fragmentation time of a jet into a virtual photon with large invariant mass M scales like 1/M. Thus 
dilepton with large mass M can be created in the medium with rather small corrections due to energy loss of their 
parent jet, while low-M dileptons should be formed outside the medium with their yield affected by the full energy 
loss suffered by the jet. Since the interesting region for the fragmentation process is for low masses M, where it is 
expected to be as important as the direct process, we assume that virtual photons fragment outside the medium after 
the parent jet has obtained its final energy. We define a medium-modified effective fragmentation function Q 



/o(z,Qf) = J d 2 r ± V(r ± ) J dE / ~-jg-(E f] E i )lJ°. /q (z',Q F ) (66) 



where z — px/Ei and z' = px/Ef. dN 3 q C q /dE is the probability to for a given quark with final energy Ef when the 
initial jet is a particle of type c and energy Ei, given by the solutions to Eq. (|50|) . ^ 7 */ 9 is the leading order vacuum 
fragmentation function from jHlai. The factor V(r ± ) is introduced to take care of the spatial distribution of jets 
and their propagation length in the medium. 
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We implement our cuts for the Drell-Yan process as well, so that the final yield is given by 



dN m _ 2 (AT, 



coll/ 



dMHy d 



OO 

/ ^^[ K ^ d^ZTd^ K, ^ d^^} p ^ - vcuupt) - (67) 



where Kdir and Kf rag are introduced to take care of higher order effects. We reproduce the Drell-Yan isolated 
muons from v / J/vtv=630 GeV p + p data at CERN 0], for low invariant masses (M < 2.5 GeV) and high transverse 
momentum (px > 6 GeV), with a constant factor Kdir — 1.5. We also take Kf rag = 1.5, quite in line with the 
A'-factor used for real photons from jet-fragmentation in Ref. Q. While we assume the AT-factors to be constant, 
they should be in principle ^/sm n, P t and M dependent. However, an estimate of the variation of Kdir with those 
parameters can be found in Ref. j3Jj. For example, next-leading order effects for the direct Drell-Yan contribution at 
y/sWJj = 2 TeV in the range 10 < p T < 20 GeV would varywithin 10% from M ~ 2.5 GeV to M ~ 4.5 GeV. 

Here we assume (N coll ) = 975, a in = 40 mb for RHIC g| and (N coll ) = 1670, a m = 72 mb for the LHC 0. In 
our calculations we use CTEQ5 parton distribution functions 0] with EKS98 shadowing corrections 0] • 

The main background at RHIC energies for the dilepton production processes considered so far is expected to 
be decay of open charm and bottom mesons 0. During the initial hard scattering, cc (bb) pairs are produced 
and can thereafter fragment into D(B) and D(B) mesons. We consider here only correlated decay, which happens 
when a positron coming from the semileptonic decay of a D(B) is measured together with the electron from the 
scmileptonic decay of a D(B). The results for heavy quark decay have been obtained with the techniques of Ref. |4f5j . 
Collisional energy loss of the heavy quarks propagating in the QGP has not been included. Since this constitutes 
a background to our process, we adopt this conservative point of view. 



VI. RESULTS 



We choose the same parameterization of the plasma phase that was previously used successfully in the studies 
of high Q and low to intermediate px photons pSlj. For central Au+Au collisions at RHIC we have an initial 
temperature Ti=370 MeV and initial time Ti=0.26 fm/c for the plasma phase, corresponding to the particle rapidity 
density dN/dy =1260 This latest value is obtained from the measured pseudo-rapidity distribution of charged 

particles dN c h/drj 0: dN/dy ~ ^\drj/dy\dN c h/dri, where \dr]/dy\ ~ 1.2 around y — 0. For LHC, our initial conditions 
are T^=845 MeV and Ti=0.088 fm/c corresponding to dN/dy=5625 F° r the processes involving the QGP, we 

assume three light flavors and we fix a s — 0.3. As the jets are defined to be particles having a transverse momentum 
greater than a scale pq, with pq 3> lGcV, we have set the dilepton momentum cutoff pr cut high enough to avoid any 
sensitivity to the choice of pq. We take pT cut = 4(8) GeV for RHIC (LHC). The cuts on leptons rapidities emulate 
the PHENIX experiment at RHIC, y C ut=0.35, while we use y cu t =0.5 at LHC. 

The dileptons produced at RHIC by the interaction of jets with the medium are shown in Fig.0 Dileptons from jet- 
pole interactions, i.e. from annihilation of a jet parton with a (g_)-mode, are negligible, while the jet-pole interactions 
involving (q + )-modes tends toward the Born term at high invariant mass, as it was the case for thermal-thermal 
reactions in Fig. |SJ On the other hand dileptons from jet-cut interactions do not behave like the cut-pole process 
in Fig. |SJ They become the dominant contribution at high- M. The expressions for the cut-pole process, shown in 
Fig. Et>+c, involve the functions /3±(|g|) with /3±(|q|) — > for \q\ 3> T. When M is large, in order to keep the value 
of \q\ modest, the energy E\ of the incoming parton has to be large: E\ > M 2 /4\q\. As the thermal phase space 
distribution function fpo decreases exponentially for large E\, the cut-pole contribution turns out to be negligible 
for large M. However, when the incoming particle is a jet with a power-law distribution, high values of E\ are 
not suppressed and the cut-pole contribution is important. Therefore the sum of all the jet-thermal processes with 
HTL effects included (solid line), is more important than the jet-thermal contribution without HTL effects (double 
dot-dashed line) by more than a factor 2 for M above 8 GeV. For M below 1 GeV HTL corrections increase the Born 
term by one order of magnitude, because of the 1/M 2 behavior in Eq. (|52|) . We have verified that when extrapolated 
to M = 0, the cut-pole contribution reproduces the result for the jet-photon conversion in the QGP Q, such that 
E1 " ^cut-pole (P. M -> °)/ d3 P = E-<dR iet _ th (p)/d 3 p. 

In Fig. El we show the mass spectrum of dileptons in central Au+Au collisions at RHIC (^/snn — 200 GeV) and 
central Pb+Pb collisions at LHC (^/snn — 5.5 TcV). For both collider energies, the jet-thermal contribution exceeds 
the thermal dilepton production by an order of magnitude, which was also the case for high-pr photon production 
However, at RHIC the dominant sources of dileptons for M > 3 GeV are heavy quark decay and the direct Drell-Yan 
process. At intermediate mass, between 1 GeV and 3 GeV, the jet-thermal contribution is as important as these 
two processes. Below 1 GeV, the fragmentation of jets into virtual photons turns out to be comparable to the direct 
Drell-Yan production and the jet-thermal contribution. At LHC, the whole range of invariant mass is dominated by 
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FIG. 7: (Color online) High-momentum dileptons produced from the interaction of jets with QGP in Au+Au collisions at 
200 GeV (RHIC). The initial temperature is Ti=370 MeV. Short-dashed line: interaction of jets with poles with positive \; 
dot-dashed line: interaction of jets with poles with negative \; long-dashed line: interaction of jets with cuts; solid line: sum 
of the latter processes; and double dot-dashed line: interaction of jet with QGP without HTL effects (Born term). See text for 
details. 
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FIG. 8: (Color online) Sources of high-pr dileptons in central Au+Au collisions at RHIC (left) and Pb+Pb at the LHC (right). 
Solid lines: semileptonic decay of heavy quarks; dashed lines: direct Drell-Yan contribution; dotted lines: jet-thermal interaction 
with HTL effects; double dot-dashed lines: jet-fragmentation process; and dot-dashed lines: thermal dilepton production with 
HTL effects. 



charm decay, but on the other hand the jet-thermal lepton pairs exceed the direct Drell-Yan yield below 7 GeV. We 
want to stress that energy loss of heavy partons has not been included here, so that the heavy quark contribution 
is likely an upper limit of what is to be expected. Lepton pair production from jet-plasma interactions being as 
important as this upper limit background at RHIC for 1 < M < 3 GeV, we thus expect the experimental detection of 
a jet-plasma contribution to be feasible. For real photons, the inclusive yield is largely dominated by the background 
coming from the decay of neutral mesons |53| . in the region where the jet-thermal contribution is expected to be 
maximal. 

While the thermally induced dilepton yield may be strongly affected by initial conditions like temperature and 
thermalization time, this is not the case for the jet-plasma contribution. This is due to the phase-space distribution 
of jets, which has a weak dependence on temperature and makes the dilepton production less pronounced in the early 
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M (GeV) 



FIG. 9: (Color online) Effect of the cut on leptons rapidities for direct Drell-Yan (top) and thermal dileptons (bottom) without 
HTL effects at RHIC. The dotted lines include all leptons from the respective process while the solid lines include only leptons 
having absolute value of rapidity smaller than 0.35. 




FIG. 10: (Color online) Effect of jet energy loss in jet-thermal production of dileptons at RHIC, with HTL effects (top) and 
without (bottom). The solid lines include the effect of energy loss, while the dashed lines don't. The results for jet-plasma 
interactions with HTL have been rescaled by a factor 10 for clarity. 



time of the QGP evolution (see Ref. |j] for a more detailed discussion of the effect of initial conditions on real photon 
production). The same argument can be also used to estimate the effects of a possible chemical non-equilibrium on the 
jet-plasma contribution. On one hand, lower quark fugacities suppress the dilepton emission at a given temperature, 
but on the other hand, smaller fugacities would imply larger temperature, thus increasing the dilepton yield. The 
interplay of those effects was studied in Ref. [54| for thermally induced dileptons, showing a suppression in the yield for 
large invariant masses at RHIC (a factor ~ 2 of suppression for M ~ 4 GeV). However, for jet-medium interactions the 
fugacities would enter into the production rate only linearly rather than quadratically for thermally induced dileptons, 
thus the effect of chemical non-equilibrium should be less pronounced. 

The effect of the momentum cut on single lepton rapidity is shown in Fig. [5] at RHIC energy for Drell-Yan and 
thermal-thermal processes without HTL corrections. For both cases, the cut reduces the yield by a factor ~ 3 and 
is almost independent of M, except in the low mass region. When M is small, the lepton rapidities tend to be very 
close to the pair rapidity yd = 0, making the cut less important. 

The effect of energy loss on the jet-thermal lepton pair production is explored in Fig. ^| for RHIC. We observe 
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that for the case without HTL effects (Born term), dileptons are reduced by about 30% for M ~ 1 GeV, while the 
suppression decreases with increasing invariant mass, reaching about 15% above M = 4 GeV. For a given invariant 
mass M and jet energy E±, the minimum energy for the thermal parton is Ei min = M 2 /AEi. This minimal value is then 
favored by the steep thermal spectrum, leading to a dependence of the yield that is roughly given by exp(— M 2 / AE\T) . 
This implies that dileptons with large mass M are more likely to be emitted at times where the temperatures T is still 
high which favors small jet propagation time and small energy loss. Including HTL effects leads to a suppression of 
about 30% for the range dominated by the cut-pole contribution, M < 2 GeV and M > 7 GeV. The region 2 < M < 7 
GeV, dominated by the pole-pole contribution, shows a weaker suppression, equivalent in magnitude to that of the 
Born term in the corresponding invariant mass range. The suppression factor is weaker for dileptons from jet-medium 
interactions than for high-p^ pions, where a large suppression factor (^4 — 5) has been observed at RHIC 0]. This 
can be understood because dileptons can be produced at any point during the propagation of the jet in the medium, 
whereas pions, due to their large formation time, are produced after a jet has left the medium, and thus they suffer 
from the full loss of energy 




FIG. 11: (Color online) pr distribution of dileptons, integrated in the range 0.5 < M < 1 GeV, for RHIC (left) and the LHC 
(right). Dotted lines: direct Drell-Yan process; double dot-dashed lines: jet-fragmentation process; solid lines: jet-thermal 
process with HTL effects; dashed lines: jet-thermal process without HTL effects; and dot-dashed lines: thermal induced 
reactions with HTL effects. 



It is also interesting to discuss the dilepton yield as a function of the dilepton transverse momentum pt in certain 
windows of the mass M. This is done by substituting the integral over pt in Eqs. (|45|l and iJBTJ, by J dM 2 / (2ttpt) ■ 




P x (GeV) p T (GeV) 

FIG. 12: (Color online) pr distribution of dileptons, integrated in the range 0.5 < M < 1 GeV (left panel) and 1.5 < M < 2.5 
GeV (right panel), for central Au+Au at RHIC. Solid line: sum of Drell-Yan process, jet-thermal and thermal-thermal process 
with HTL effects; dashed line: sum of Drell-Yan process, jet-thermal without HTL and thermal-thermal with HTL effects; and 
dot-dashed line: sum of Drell-Yan process and thermal-thermal reactions with HTL effects. 
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The results for RHIC and the LHC are displayed in Fig. ^2 f° r the mass integrated in the range 0.5 < M <1 
GeV. The ordering of the contributing sources here is very similar to the one seen for real photons as a function of 
Pt in Ref. The direct component of the Drell-Yan process dominates for high-py dileptons at RHIC while the 
jet-thermal contribution, with HTL, dominates for px <5 GeV. The hard thermal loop calculation enhances this yield 
by more than a factor 4 compared to the jet-thermal contribution without HTL. At the LHC, jet-thermal dileptons 
(HTL effects included) are the most important source in the entire px range, 8 < px < 17 GeV. The jet-thermal 
interaction appears to be as important for dileptons as it was for real photon production. 

The total direct dilepton spectrum for RHIC is shown in the left panel of Fig. The solid line includes Drell- 
Yan and QGP contribution (jet-thermal and thermal-thermal) with HTL effects. Leaving out the HTL resummation 
for jet-thermal dileptons (dashed line) reduces the yield by a factor ~ 1.5 around pt=4 GeV. The absence of any 
jet-thermal interactions at all (dot-dashed line) would reduce the total yield by a factor ~ 2 at px—^ GeV. This 
emphasizes the importance of this process in the presence of a QGP. 

The only potentially important contribution that is not included in our work is in-medium bremsstrahlung (qi — > 
qij*) and annihilation (qqi — ► of an incoming thermal parton or jet, where i denotes a quark, antiquark or 
gluon. This goes beyond the current formulation of AMY. However, they have been calculated in Ref. |5J| for the 
case of incoming thermal partons, showing that for low mass dileptons, the bremsstrahlung and annihilation more 
than double the contribution obtained from 2 — > 2 processes (q + g — > q + 7*, q + q — > g + 7*). It thus turns out 
that thermally induced bremsstrahlung and annihilation are as important relatively to 2 — > 2 processes, no matter if 
the photons are virtual or real |56|. On the other hand, for the case of real photons and incoming jets, it has been 
shown in Ref. 4} that those bremsstrahlung and annihilation processes are reduced by a factor 3-4, relatively to the 
in-medium jet-photon conversion process. This is because the momentum distribution of jets is less steep than the 
thermal one, making the jet-photon conversion the most important in-medium process. Therefore, if one assumes 
that the in-medium jet-bremsstrahlung contributes in the same way to virtual and to real photons, it would enhance 
the solid line by less than 15%. 

Finally, the right panel of Fig. H2l shows the dilepton spectrum for an another invariant mass window. It is located 
at higher values 1.5 < M < 2.5 GeV between the (j> and the J/ip masses. As we could have expected from Fig. 
the effect of the HTL resummation is not very important for this mass window. However interactions of jets with the 
plasma are still a very important source of dileptons, and this should be detectable. 

VII. SUMMARY AND CONCLUSIONS 

In this work we presented calculations of different sources of lepton pairs in high energy nuclear collisions. We 
take into account Drell-Yan, fragmentation from jets, QGP contributions and heavy quark decay. Hard thermal loop 
resummation has been included in the calculation of the leading order photon self-energy. We explicitly checked 
that the imaginary part of this self-energy, evaluated within finite-temperature field theory, and a different approach 
starting from relativistic kinetic theory and using the corresponding Feynman amplitudes, lead to the same results. 
We obtain the jet-plasma interaction by substituting the phase-space distribution of one incoming thermal parton, 
by the distribution of jets. While the HTL corrections are important for thermal-thermal processes at low-invariant 
mass, they are important for both low and high invariant mass when the incoming parton is a jet. 

For the high momentum window investigated here, px > 4 GeV, dilepton emissions due to jet-plasma interactions 
are found to be much larger than thermal dilepton emission. At low to intermediate dilepton mass, productions 
from jet-plasma interactions are comparable in size to the Drell-Yan contribution and constitute a good signature for 
the presence of a quark gluon plasma provided the dominant background of heavy quark decay could be subtracted. 
The AMY formalism has been used to account for energy loss of jets in the QCD plasma; this energy degradation 
reduces the effect of jet-thermal processes by ~ 30%. Further study involving heavy quark energy loss will be needed 
to obtain a better estimate of this channel, together with an explicit calculation of dileptons from medium-induced 
bremsstrahlung. 
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